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We present a model for the disorder of the heterointerfaces in GaAs quantum wells including 
long-range components like monolayer island formation induced by the surface diffusion during the 
epitaxial growth process. Taking into account both interfaces, a disorder potential for the exciton 
motion in the quantum well plane is derived. The excitonic optical properties are calculated using 
either a time-propagation of the excitonic polarization with a phenomenological dephasing, or a 
full exciton eigenstate model including microscopic radiative decay and phonon scattering rates. 
While the results of the two methods are generally similar, the eigenstate model does predict a 
distribution of dephasing rates and a somewhat modified spectral response. Comparing the results 
with measured absorption and resonant Rayleigh scattering in GaAs/AlAs quantum wells subjected 
to growth interrupts, their specific disorder parameters like correlation lengths and interface flatness 
are determined. We find that the long-range disorder in the two heterointerfaces is highly correlated, 
having rather similar average in-plane correlation lengths of about 60 and 90 nm. The distribution of 
dephasing rates observed in the experiment is in agreement with the results of the eigenstate model. 
Finally, we simulate highly spatially resolved optical experiments resolving individual exciton states 
in the deduced interface structure. 
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I. INTRODUCTION 



The interfaces of a semiconductor quantum well (QW) 
always present some amount of disorder as a natural con- 
sequence of the physics of the growth process. This dis- 
order determines to a large extent the excitonic optical 
response both in frequency and time domain. Optical 
spectroscopy on QW excitons is therefore able to probe 
the interface disorder. However, the link between the 
optical response, particularly resonant Rayleigh scatter- 
ing (RRS), and the details of the interface structure, is 
not straightforward and calls for a microscopic quantum 
model. QW interface fluctuations affect both the po- 
tential and the kinetic contributions to the total exciton 
center-of-mass (COM) energy, resulting in a Schrodinger 
equation for the exciton COM wavcfunction along the 
QW plane, including a static disorder potential that is re- 
lated to the actual interface structure. The exciton COM 
wavcfunction is proportional to the exciton macroscopic 
interband polarization and relates therefore directly to 
the linear optical response of the system. The eigen- 
states in this model present both a spatial localization 
within the QW plane and an inhomogcneous distribu- 
tion of cigen-cnergics. When QWs arc grown by molec- 
ular beam epitaxy, interface disorder is mostly deter- 
mined by surface diffusion and segregation. When using 
semiconductor alloys for the well or the barrier, also al- 
loy disorder becomes important. While surface diffusion 
tends to create laterally extended regions with a well- 
defined monolayer (ML) thickness, segregation and al- 
loy disorder result in nanoroughness on an atomic length 
scale, smaller than the exciton Bohr radiusi. The in- 



terruption of the epitaxial growth at the QW interfaces 
leads to a longer surface diffusion length, creating ML 
islands in the 100 nm size range that have been observed 
in various experimental studiesSiii. When ML islands 
are sufficientljii large compared to the excitonic Bohr ra- 
dius, the exciton spectrum can split into several peaks 
of reduced inhomogcneous broadening^. In a simple pic- 
ture these peaks are resulting from exciton states present 
within flat regions of different integer ML thickness, while 
the inhomogcneous broadening results from a residual 
nanoroughness. To create from this picture a realistic 
description of the ensemble properties of the exciton in 
the QW, one needs to consider the following issues: (i) 
the statistical distribution of size and shape of regions 
with different ML thickness of the interfaces; (ii) the dif- 
ference of the structural parameters and the statistical 
correlation of the two interfaces; (iii) the lateral quan- 
tum confinement energy of the exciton states; (iv) the 
coupling of the exciton states to light. Once the statisti- 
cal properties of the interface structure arc described by a 
model with a small number of parameters, these param- 
eters can be deduced from the measured RRS spectrum 
and dynamics. 

Existing theoretical models of RRS in QWs were of- 
ten based on very simplified assumptions for the shape 
of the disorder potential. The simplest picture, intro- 
duced in the original work on RRS theory by Savona and 
Zimmermanni consisted in assuming a random poten- 
tial with Gaussian or exponential correlation length in 
space. Although this approach can explain the principial 
properties of RRSf^ it is unable to predict the dynam- 
ics of more complex situations like the one considered 
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here. In a previous workS, we had developed a rudimen- 
tary extension of the disorder model including a short- 
and a long-range disorder component, for the modeling 
of RRS in samples where one of the two heterointerfaces 
was growth interrupted. 

In this work we apply the microscopic theory of the 
RRS optical response of cxcitonsi to a situation where a 
pronounced ML splitting is present. To this purpose, we 
develop a comprehensive model of the disorder potential 
based on a separate description of the direct and inverse 
heterointerfaces of the QW. The model accounts for a 
ML island structure at both QW interfaces, with differ- 
ent correlation lengths, anisotropy factors and statistical 
distributions of the integer values of the ML thickness. 
It also includes nanoroughness and a parameter that ex- 
presses the statistical correlation between the two inter- 
faces. The experimental RRS data, both spectrally and 
time-resolved, were taken on a GaAs/AlAs QW sample 
grown with long growth interruptions at both interfaces. 
The amount of deposited GaAs is varying slowly across 
the lateral position on the sample, so that a fine control 
of the average QW thickness by a small fraction of one 
ML is achieved. This in turn produces a continuous evo- 
lution of the ML-split exciton spectrum, which replicates 
itself after a thickness change of an integer ML with only 
an energy shift given by the corresponding difference in 
QW confinement energy. Correspondingly, also the time- 
resolved RRS shows a peculiar dependence on the frac- 
tional ML thickness. We use this rich phenomenology to 
determine the parameters of a detailed interface model 
that we develop. The calculated RRS response is con- 
fronted with the measured spectral and time-dependent 
RRS signatures, and the effects of the various interface 
properties are elucidated. This analysis allows an estima- 
tion of the actual parameters characterizing the investi- 
gated sample - to a much larger extent than in previous 
studies based on the same RRS theory. Specifically we 
find an almost perfect correlation between the direct and 
inverse QW interfaces ~- a novel result that is of general 
relevance for epitaxial QW samples. In the second part 
of the work, we turn to the numerical computation of the 
actual exciton COM eigenstates as resulting from the dis- 
order model. This part of our study is required in order 
to account for state-dependent radiative recombination 
and phonon scattering rates, which prove important es- 
pecially in the RRS dynamics. It provides us with further 
insights into the microscopic structure of the single exci- 
ton states and their optical properties. 

The paper is organized as follows. In Section II we de- 
scribe the time-propagation RRS theory and the interface 
disorder model, followed by a study of the dependence of 
the simulated RRS on the model parameters. Section III 
describes the experimental data and the features of the 
QW interfaces, as deduced from the comparison to the 
model calculations. Section IV is devoted to a micro- 
scopic simulation of the exciton eigenstates and of their 
radiative and phonon scattering rates. Results for en- 
semble averages and individual disorder realizations are 



presented. Section V contains the conclusions and out- 
look of the work. 



II. THEORETICAL DESCRIPTION 

The present model is based on the theory of RRS 
within the exciton COM approximation, as derived by 
Savona and Zimmermanni. The exciton interband polar- 
ization P is proportional to the exciton COM wavefunc- 
tion along the QW plane and therefore evolves according 
to the time-dependent Schrodinger equation 

-f ^^in(R,t), (1) 

where R = {x, y) is the position in the QW plane, M 
is the in-plane exciton mass, /x is the dipole matrix el- 
ement of the excitonic transition, and Ei^ (R, t) is the 
resonant excitation field. In the present model, we ne- 
glect the vector character of the electric field and of the 
interband polarization. Rayleigh scattering is expected 
to conserve the orientation of the polarization vector of 
the scattered field, whereas the long-range exchange part 
of the Coulomb interaction can induce a rotation of this 
vector on a timescale of several tens of picosecondsiSiii. 
The static COM potential T^(R) is an effective quantity, 
related to the actual exciton energy fluctuations via a 
convolution with the Is-exciton relative wavefunctionii. 
Differently from the usual deflnition, we have incorpo- 
rated here the average Is exciton energy Ex into 1^(R). 
We detect the RRS in the far field and neglect time delays 
due to light-propagation. In linear response theory, the 
RRS amplitude i?out (k, t) is proportional to the Green's 
fimction G(k, kin,<) of the Schrodinger equation ex- 
pressed in momentum space^iii where kin is the mo- 
mentum of the incident electromagnetic field ii^in(kin, t), 
assumed to be a plane wave. To determine the spectrally 
resolved RRS Kut(k,w) = Kut(k, t)e(^'^-^)*dt, we 
introduce a phcnomcnological polarization decay rate 7. 
The resulting RRS intensity 

/RR,s(k,t) cx \Eout{Kt)\'^ , ^RRs(k,a;) cx |£'out(k, a;)|^ 

(2) 

is averaged over a large statistical ensemble of disorder 
realizations. In a similar way, we can define the exciton 
COM spectral function A(k,w) = -2Im{G(k, k, w)} as 
proportional to the imaginary part of the COM Green's 
function at equal momentum arguments. The spectral 
function at momentum kin, averaged over disorder real- 
izations, is proportional to the exciton absorption spec- 
trum (the same average is identically zero if the Green's 
function is taken with two differing momentum argu- 
ments). The simulated spectral function will be com- 
pared, in this and in the next Section, to the measured 
excitonic absorption. Throughout the paper, when quot- 
ing values for frequencies or rates in energy units, they 
arc multiplied by h. 
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A. Model for the interface structure 

In order to model the exciton COM potential in pres- 
ence of ML islands in the heterointerfaces, as in the QWs 
studied here, the following procedure is used. First, a 
random white- noise function Wq(R) is generated, such 
that (u;a(R)wa(R')> = (5(R- R'). The index a = 1,2 
denotes the two QW heterointerfaces. We then derive a 
spatially correlated function of unitary standard devia- 
tion and zero mean value 



Wa (R) (X / dR' Wa (R') exp 



{x-x'f {y-y'f 



(3) 

Here we have introduced the Gaussian correlation lengths 
with J = x,y, which will determine the length scale 
of the ML islands along the x and y directions, respec- 
tively. The separate length scales for the two directions 
allow to describe a spatial anisotropy of the ML struc- 
ture. We assume equal orientation of the anisotropy axis 
for both interfaces since we expect the microscopic ori- 
gin of the anisotropy to be given by the anisotropic re- 
construction of the [001] growth surface along the [110] 
in-plane direction, equal for the whole epitaxially grown 
structure. We express the two correlation lengths as 
a function of an anisotropy parameter and of an 

1 /2 

isotropic correlation length ^q,, defining £^x,a = Ca £,a 

— 1/2 

and ^y^a = €a ^a- We proceed by generating a ML step 
function from Wa that models the property of a growth 
surface to form compact regions of integer ML coverage 
due to the lateral binding energyi^ of the adatoms, as 
opposed to a random arrangement. We use the function 



la (R) = Round 



Waj-R) 

Co 



(4) 



where the function Round [.t] denotes the nearest integer 
of a real number x. The quantity C.a is the interface flat- 
ness, controlling the height variance of the interface to be 
approximately C,~^ MLs. The parameter A' e [0, 1] is a 
fractional ML shift that describes the fractional ML cov- 
erage of the first interface, i.e. the AlAs surface on which 
the GaAs QW is grown. The spatial average of the ML 
coverage (/i)r is approximately A' for interfaces that 
consist of several ML steps, i.e. C ^ 3. This is the case 
for the simulations presented later on. For flatter inter- 
faces the fractional ML coverage is a non-linear function 
of A' (see upper inset of Fig.QJ. The same arguments 
hold for (/2)r and A' + A. Due to an unavoidable slight 
misorientation (>1") of the substrate with respect to the 
[001] direction, the fractional ML coverage of the inter- 
face is expected to vary over many ML thicknesses across 
the excited surface area, which has a lateral extension of 
about 10^ ML (300 iim) . We will therefore use a uniform 
distribution of A' over [0, 1] in the statistical ensembles 
for the simulations. The quantity A is shifting the second 
hctcrointcrface with respect to the first one, as specified 
by the Kronecker 5a^2 in Eq. It defines the average 



distance between the first and the second heterointcrface, 
i.e. the average QW thickness. This thickness can be ac- 
curately varied in the experiment by displacing the laser 
spot along the QW plane since the average QW thickness 
is given by the deposited amount of GaAs within the sur- 
face diffusion length and varies only by about 0.1 ML over 
300 /^m lateral extensioni^. In this argument, we assume 
that Ga and Al atoms do not evaporate from the sur- 
face once deposited. This assumption is justified for the 
growth parameters used for the investigated samples. 

We continue by deriving a local exciton energy 
[/ (R) = F (/2 (R) - h (R)) from the local ML-thickness 
/2 — /i . The dependence of the exciton energy on the ML- 
thickncss has the form 

F(A) =ao + ai/A + a^/A^ + ag/A^ + a^/A^ (5) 

where the parameters a„ are fitted to the experimen- 
tal datai^ obtained on the same sample by measur- 
ing the exciton peak position at equivalent fractional 
ML thicknesses over a large range of thicknesses (sec 
Fig.^l. A random, spatially uncorrelated potential u (R) 
is summed in order to account for short-range disorder 
originating from segregation. The COM potential V^(R) 
is obtained by convoluting [/(R) + m(R) with the elec- 
tron and hole probability distribution in the Is-excitonS 
|0(R')|^ cx exp(— |R — R'I/ob), using an exciton Bohr 
radius ob = 8 nm. This value was deduced from the 
measured diamagnetic shift ad — 33/xcV/K and the re- 
duced mass /.i* = 0.0586 for a 37ML thick GaAs/ AlAs 
quantum welli4 using the formula = ag(3e^)/(16/i*) 
valid for an in-plane exciton wavefunction proportional to 
exp(— r/ae), with the in-plane electron-hole distance r. 
The amplitude of the short-range disorder w(R) is param- 
eterized by cr, expressing its standard deviation after the 
convolution with the Is-cxciton relative wavefunction. 

Within this model, particular attention must be payed 
to the statistical correlation between the ML fluctuations 
of the two interfaces. In the epitaxial growth, a sub- 
stantial amount of correlation can be expected between 
the ML fluctuations of the lower and upper interfacesi^. 
Wc introduce this correlation into our model via the two 
white- noise functions u'q.(R) that enter Eq. (j2J). The 
limiting case of fully correlated heterointerfaces is ob- 
tained by taking w\ = wi. Conversely, the fully uncorre- 
lated case results from using two uncorrelated white noise 
functions w\ and wi. In order to model the intermedi- 
ate case of partial correlation, we introduce two mutu- 
ally uncorrelated white-noise functions u>(R) and u)'(R). 
We could then proceed by simply defining w\ — w and 
u'2 = t^w + (1 — K^^^l'^w' , therefore assuming a uniform 
correlation k for all the Fourier components of w\_i. How- 
ever, the total number of Ga atoms deposited over length 
scales larger than the surface diffusion length has to be 
conserved in the growth process. Even though the atoms 
are deposited at random in molecular beam epitaxy, the 
total number of Ga atoms contained in the QW within 
an area given by the surface diffusion length is > 10^, 
so that the statistical error in the QW thickness is be- 
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QW thickness (ML) 

FIG. 1: Measured exciton transition energy (circles) versus 
nominal QW thickness A, and the fitted F{A) (solid line) 
using (ao = 1.516, ai — —0.233, a2 — 92.5, a3 — —806, 04 — 
2739) eV. The measured exciton PL peak energy at a thickness 
resulting in a dominating single ML peak (see lower inset) is 
attributed to the exciton potential F{A) at +0.15 fi^actional 
ML thickness. This slight shift is needed to account for the 
quantum confinement due to short-range disorder, which is 
shifting the peak in the optical spectrum to slightly lower en- 
ergy with respect to the nominal zero-point (spatial average) 
of the disorder potential (see e.g. Fig. 5 in Ref. ITU) . The 
value -1-0.15 ML is giving the best fit between RRS measure- 
ments and simulations at A « 38. Additionally, the mea- 
sured ML splitting (squares) and its fit dF/dA (dashed line) 
is shown. The upper inset shows the average ML coverage of 
the first interface {Ii)r as function of A'. 



low 10^'^, a fluctuation that can be neglected compared 
to the interface roughness. The long-wavelength compo- 
nents of the ML fluctuations in the two interfaces have 
therefore to be fully correlated. In order to accommo- 
date this requirement, we modify the above approach for 
the long- wavelength component of W2 ■ Using the Fourier 
transforms and w'^ of the two generating functions, 
we define the function W2 (R) as the inverse Fourier trans- 
form of 



W2,k = fkWk + {1 - fi 



1/2. 



with 



fk = (1 - K)exp - 



(6) 



(7) 



where Ic is a surface diffusion length. In this way, the pa- 
rameter K €E [0, 1] still represents the correlation between 
the two interfaces for the Fourier components having 
wavelengths shorter than Ic, while for longer wavelengths 
the interfaces are fully correlated as dictated by the con- 
servation of Ga atoms. It the following simulations, Ic is 
always much larger than the correlation lengths ^q, and 
has therefore no significant effect. We used Ic = 500 nm. 
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FIG. 2: (a) Color plot of the first interface profile /i(R) for 
A = 38.0 while A' is varied in the interval [0, 1] along the 
y-position. The color steps represent ML height steps, (b) 
Same for the second interface profile J2(R-). (c) COM po- 
tential V'(R) for an integer monolayer thickness A = 38.0. 
(d) Same for half-integer monolayer thickness A = 38.5. The 
parameters for this potential are given in Table|l] and the 
correlation is k = 1. The color scale in c) & d) covers 5meV. 



The conservation of Ga atoms also requires that the 
long wavelength components of W1.2 created from uii.2 
are equal, which relates the correlation length and the 
flatness parameters of both interfaces by ^i/Ci = £,2/(2- 
We will assume this relation to hold throughout this 
work. 



B. Simulations 

The time-dependent RRS field is simulated by solving 
the time-dependent Schrodingcr equation with the poten- 
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tial 1^(R), according to Eq. The Fourier transform 
of the interband polarization then provides the RRS ani- 
phtude in wavevector and frequency domain^. The RRS 
intensity computed in this way is averaged over 10'^ — lO''^ 
reahzations of V^(R) to create a statistically relevant en- 
semble average, converging to its asymptotic value within 
3—1% accuracy. The remaining statistical fluctuations 
are only appreciable in the time-resolved RRS intensity 
at long times, at which the asymptotic RRS dynamics 
is very slow. It was showni that the realization average 
coincides with the speckle average which is performed 
in the experiment by using a large excitation region of 
about (200 /J.m)^ and averaging over a small span of scat- 
tering directionspi^ This correspondence is very useful, 
as simulating directly a large area is numerically pro- 
hibitive. We use typically a (500 nm)^ simulation area 
over 64 x 64 grid points with periodic boundary condi- 
tions. We assume a plane- wave excitation field at k = 
and take the scattered RRS field at the smallest available 
wavevector value k=47r//xm on the simulation grid. The 
average ML thickness A is known for each position on the 
sample by a previous studyi^. The free parameters used 
in this model are therefore the correlation lengths £,a, the 
flatness ^i, the anisotropy constants Cq., the correlation 
K, and the short-range disorder amplitude a. 

In this Section we proceed by discussing the simu- 
lated, spectrally or time-resolved RRS intensity for the 
disorder parameters that best describe our experimental 
data, which are presented in full detail in the next Sec- 
tion (see Fig.inj. These parameters were obtained by a 
trial-and-error procedure that consisted in (i) performing 
a simulation with a given parameter set, (ii) comparing 
to the measured data (both spectral and time-resolved), 
and (lit) adjusting heuristically the parameters and go- 
ing back to step (i). A standard best-fit procedure is 
made prohibitive by the long computing time required 
for a single simulation (a few hours for TP, a few days 
for ES). Starting from this parameter set, indicated in 
TableQlTP) we discuss how varying the different param- 
eters of the model affects the outcome of the simulations, 
also in light of previous reports concerning the interface 
structure in epitaxial crystal growth. This analysis also 
provides estimations of the relative accuracy of the pa- 
rameter values determined from the comparison to the 
experiments. 

Fig.EJa) and (b) show examples of the two interface 
functions /i(R) and /2(R) for the interface parameters 
found to reproduce the measured RRS (see Table TP). 
The fractional ML shift A' was varied linearly from 
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87 nm 


1.2 


1.6 


1.0 
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ES 


55 nm 


82 nm 


1.46 


2.2 


1.0 


1.0 


0.8 meV 





TABLE I: Parameters of the interface model obtained by 
adjusting the simulations to the measured RRS data for the 
time-propagation (TP) and the eigenstate (ES) calculations. 
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FIG. 3; Calculated spectral disorder potential density 
Dv{ijj), absorption (abs), and RRS intensity (RRS) for A = 
38.0 and 38.4 ML, other parameters as in Fig.|5| Data are 
vertically offset for clarity. The dotted lines represent the 
potential F{A) for the labelled integer A. 

to 1 along the y direction. These false color pictures il- 
lustrate how the next ML level gradually develops when 
moving along the y direction. The two interfaces are 
fully correlated {k ~ 1), as is reflected in the similar 
shapes of the two ML patterns. The non-equal correla- 
tion lengths ^1 < ^2 result in smaller features and more 
irregular patterns being present on the first interface (a). 
The resulting COM potential T^(R) for integer and half- 
integer ML thickness A of the QW are displayed in (c) 
and (d). Both pictures show that the variation of the 
fractional ML shift A' is compensated when subtract- 
ing the ML profiles (Eq.QJl of the two interfaces, and 
does thus not significantly affect the COM potential. On 
the other hand, the difference between integer and half- 
integer ML thickness appears clearly in the plots. Be- 
cause of interface correlation, for integer ML thickness 
the ML steps on the two interfaces appear synchronously 
and, in the limiting case of identical interfaces, the QW 
thickness would remain constant S The fluctuations in 
the QW thickness are thus due to the different charac- 
teristic sizes of ML islands on the two interfaces which 
give rise to ring and stripe patterns in the QW thickness 
and thus in 1^(R), as can be seen in Fig.|2Ic). In the 
case of half-integer ML thickness the steps on the two 
interfaces are out of phase, and the ML steps are likely 
to appear on one interface only, giving rise to flat regions 
in V^(R) as displayed in Fig.|2Id), which are representing 
two rather well defined ML thicknesses. 
The spectral density 

Dviuj) = A-^ J S{V{R.) - huj)dR. (8) 

of the resulting disorder potential is given in Fig.|3| For 
integer ML thickness A = 38.0, Dy has a dominant peak 
around F(38), with tails towards the high and low en- 
ergy extending to the potentials of the adjacent integer 
ML thicknesses. Due to the averaging over the exciton 
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FIG. 4: Calculated time- (a) or spectrally (b) resolved 
RRS intensity versus the fractional ML thickness A for a 
GaAs/AlAs interface of long and a AlAs/GaAs interface of 
short correlation length (^i = 25 nm, ^2 = 200 nm), and a 
ML flatness ("1 = 2.5. Data are offset for clarity. 



size, the regions of adjacent ML thickness do not result 
in such distinct peaks since they have a small spatial 
extension. The calculated absorption peak is shifted to 
higher energies due to the finite mass of the exciton giv- 
ing a lateral quantization energy. Furthermore, only the 
lower energy ML is visible in the absorption, but not the 
higher one. This asymmetry relates to the small momen- 
tum of the light, which essentially probes the fc = com- 
ponents of the wavcfunctions. With increasing energy, 
the wavefunctions get more nodes, i.e. a smaller k = 
component, so that only the lower part of the potential 
landscape relates to optically active exciton states. The 
RRS spectrum shows a similar behaviour as the absorp- 
tion, but in a more pronounced way. This is expected 
because the RRS intensity - being proportional to the 
squared COM Green's function - varies as the fourth 
power of the COM wavefunction at small momentum, 
as opposed to the second power dependency characteris- 
ing the absorption spectrum. For about half-integer ML 
thickness A — 38.4, Dv{oj) is distributed between F(38) 
and F{39), as expected. The absorption and the RRS 
shows two peaks, with the lower energy peak dominat- 
ing, again due to the dominant fc = component of the 
lowest energy eigenstates. 



1. Long correlation length of GaAs/AlAs interface 

The parameters of Table define fully correlated in- 
terfaces with similar correlation lengths ^1 — 65 nm and 
^2 = 87 nm. From structural investigations, it could be 
expected that the islands on the GaAs surface are much 
larger than on the AlAs surface due to the largely differ- 
ent surface mobility at the growth temperature of 630° C. 
On GaAs surfaces after 60 s growth interruption, islands 
of several hundred nanometers were observed^*^. Earlier 
investigations, however, measured the GaAs surface after 
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FIG. 5: Calculated time- (a) or spectrally (b) resolved RRS 
intensity versus the interface correlation k. Other parameters 
as in tabled] Data are offset for clarity. 



cooling, which can have a significantly larger correlation 
length than interior intcrfacesii. In order to test this ex- 
pectation, we therefore discuss now the RRS spectra of 
excitons in interface structures with a large correlation 
length of the GaAs to AlAs interface of ^2 = 200 nm. 
This is much larger than the localization length of a ML 
island^ Co = nh/y/2M{dF/dA) ss 18 nm for A = 38. 
This interface thus does not significantly influence the ex- 
citon localization. The free parameters are then the AlAs 
to GaAs interface correlation length d and flatness Ci. 
The variation of the RRS spectra with increasing A for 38 
to 39 is best reproduced using d — 25 nm and = 2.5, 
as displayed in Fig. 01 The resulting RRS dynamics, how- 
ever, is essentially independent of A, in contrast to the 
experimental results (see Fig.^. This finding is insensi- 
tive to the specific values of d and Ci, and shows that 
two interfaces with largely different correlation lengths 
are not consistent with the measured RRS dynamics. 
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2. Interface Correlation 

For the case of largely different lengths ^i_2 of the two 
interfaces as discussed in the previous Subsection, the 
correlation k, is not important since the ML islands of 
the interfaces arc very different. For similar correlation 
lengths of the two interfaces instead, k is important, since 
for full correlation (k = 1) only the small difference be- 
tween the two correlated interfaces remains, leaving only 
narrow flat regions in V(R) as seen in Fig. EJc). With- 
out correlation (k = 0), on the other hand, the difference 
between the uncorrelated interfaces is similar to the case 
of a single interface with a reduced flatness. Starting 
from the parameters in Table HfXP), we show in Fig.^l 
the effect of a decreasing correlation on the RRS spec- 
tra and dynamics. Already for a small deviation from 
full correlation [n = 0.95), for which the RRS spectrum 
is not signiflcantly changed, the RRS dynamics devel- 
ops an initial peak indicating the presence of a disorder 
component of large correlation lengtl^i. This component 
results from the difference between the k > \ Fourier 
components of Wi and W2, which does not vanish for 
non-perfect correlation k < 1. For even smaller values 
of K, the RRS spectrum gradually develops more ML 
peaks (A = 37,39), and the initial peak in the RRS 
dynamics gets more pronounced. The RRS dynamics 
becomes essentially independent of the fractional A, a 
behavior inconsistent with the experiment. We thus con- 
clude that the ML steps on the two interfaces have to 
be nearly perfectly correlated. This finding suggests that 
for our samples the surface diffusion length Ic is actually 
similar to the correlation lengths ^q,, so that the almost 
perfect interface correlation is created by the conserva- 
tion of the deposited Ga atoms. We would expect this 
finding to depend on the duration of the growth inter- 
ruption at the GaAs to AlAs interface and the growth 
temperature. For long growth interruption times and in- 
creasing temperature, the surface diffusion length is in- 
creasing compared to the island size, as the surface dif- 
fusion distance is increasing with \/t, while the island 
size is thermodynamically limited. With increasing du- 
ration of the growth interruption (GI) and growth tem- 
perature, the interface correlation is therefore expected 
to be diminished. This is actually supported by previ- 
ous experiments. Yu et ali^^ have studied the PL spectra 
of GaAs/ AlAs QWs depending on the GI time. They 
observe (Fig. 1 of Ref. [H) the formation of the ML 
splitting during the first 60s GI, with only 2 ML peaks 
present, while for 120 s GI, already 4 peaks are visible, 
similar to our simulation for k = 0.5 in Fig.|Sl For 
GI GaAs QWs with mixed crystal barriersiSiSSiSi nor- 
mally more than 2 ML peaks arc present. This indicates 
that for mixed crystal barriers, the interface disorder is 
partly uncorrelated, as was also observed in structural 
investigations^^. This finding can be understood consid- 
ering that a Ala;Gai_a;As mixed crystal growth surface 
consists essentially only of Ga atoms that are mobile. 
The Al distribution, which determines the interface cre- 
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FIG. 6: Calculated time- (a) or spectrally (b) resolved RRS 
intensity versus the flatness ("2 for C1/C2 = 0.75. Results for 
A = 38.0 (left) and A = 38.4 (right) are shown. Data are 
ofi'set for clarity. 

ated after overgrowth with GaAs, is therefore essentially 
unmodified during the growth interrupt. The resulting 
first interface is consequently dominated by short range 
disorder, and the ML islands formed by the growth inter- 
rupt at the second interface give rise to a larger number 
of ML peaks due to the missing correlated islands on the 
first interface. 



3. Influence of flatness and correlation lengths 

We now discuss the influence of the flatness (i and the 
correlation lengths ^1.2 on the optical response, starting 
from the parameters in Table TP). The correlation is 
taken to be unity (k = 1), and C2/C1 = ^2/£,i- We first 
change interface flatness, from the value Ci = 1-6, to a 
smoother Ci = 2.5, and rougher (Ci = 1.0) interface. The 
calculated RRS is shown in Fig.EI Generally, increas- 
ing roughness leads to a spectral broadening of the ML 
peaks, while their number is hardly affected. In the time- 
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domain, increasing roughness leads to a decrease of the 
RRS rise time, and to a suppression of the beating be- 
tween the ML peaks. From the obtained dependency, we 
estimate the accuracy of the value in TablcljTP) for the 
flatness Ci to about ±30%. 

Let's now look at the influence of the correlation length 
^1 while keeping the ratio = 0.75 constant. The 

corresponding simulations are displayed in Fig.[7| In- 
creasing this length decreases the lateral quantization 
energy in the ML islands. At ^1 = 35 nm, the ML 
peaks are not developed, since the quantization energy 
is comparable to the ML splitting. With increasing ^1, 
the ML peaks develop into well defined resonances. In 
the RRS dynamics, this trend is reflected in the appear- 
ance of a well defined temporal beating between the ML 
peaks. Note the absence of an initial RRS peak pre- 
dicted for longer correlation length in a simple disorder 
potentiaii^. This is related to the perfect interface cor- 
relation (compare also with Fig.|3J), effectively canceling 
the long-wavelength components of the disorder. From 
the obtained dependency, we can estimate the accuracy 
of the value of the correlation length ^1 in Table QlTP) 
to about ±15%. 



photon energy (eV) 

FIG. 8: Calculated time- (a) or spectrally (b) resolved RRS 
intensity versus the correlation length ratio ^1/^2 for = 
75 nm. Results for A = 38.0 (left) and A = 38.4 (right) are 
shown. Data are vertically offset for clarity. The dashed and 
dotted lines in the right lower panel at ^1/^2 = 1 are for 
A = 38.2 and A = 38.6. 



Finally, in Fig. |H1 we investigate the influence of the 
ratio of the correlation lengths ^1/^2 for a fixed value of 
^1*^2- For unity ratio, a difference of the interface struc- 
ture exists in our model only for non-integer ML thick- 
ness. For integer A the interface difference is constant, 
resulting in a spectrally sharp RRS only broadened by 
the short-range disorder a. Varying A continuously from 
one integer value to the next, rings and stripes of larger 
thickness form, gradually getting wider until they en- 
tirely cover the plane as A reaches its next integer value. 
The rather well defined ring width leads to a third peak 
in the spectrum, with an energy above the 39 ML peak 
given by the one-dimensional quantization energy across 
the ring. The peak accordingly shifts to lower energy 
with increasing A. 

With increasing deviation of ^1/^2 from unity, the ef- 
fect of the interface correlation decreases, so that for 
^1/^2 = 0.5 more than two ML peaks are visible, and 
the RRS dynamics shows an initial peak, similar to the 
results for decreasing correlation in Fig.|Sl and develop- 
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ing into the case shown in Fig.0J From the obtained 
dependency, we can estimate the accuracy of the value 
in Tablc|Il;TP)for the ratio ^1/^2 to about ±8%. 



4-. Influence of anisotropy and short-range disorder 

The influence of the correlation length anisotropics 
within 0.5 < eo, < 2 on the simulation results is not 
very pronounced, and is similar to the combined effects of 
other parameters. We thus have chosen the isotropic case 
= f in the following. Wc remark that an anisotropy 
is expected due to the surface reconstruction. The effect 
of such an anisotropy on the exciton finc-structure was 
observedS^, and an analysis in terms of a correlated dis- 
order potential in a 5nm wide GaAs/Alo.sGao.yAs QW 
grown without GI revealediS a value of e « 0.6. 

The short range disorder described by a leads to an 
overall spectral broadening. Its main effect for the pa- 
rameter set in Tabic UTP) is to broaden the sharp ML 
peaks developing for A close to an integer ML thickness. 
The present value is 0.8 meV. To interpret this value, we 
consider that the volume over which the exciton averages 
contains about N = 30000 Ga/Al sites. For a random 
distribution of n Al sites {n <^ N), we would expecli^ 
a disorder potential variance of cr = ^/nEc/N, where 
Ec = 1.6 eV is the band gap difference between GaAs 
and AlAs at the F point. From a = 0.8 meV we estimate 
an average Al concentration of n/N « 0.75%. This rough 
estimate is consistent with the expected segregation^^. 



III. COMPARISON WITH THE EXPERIMENT 

We now compare the calculated optical properties with 
the measured RRS and optical density of the sample de- 
scribed as wafer 5 in Ref.|3 K contains single GaAs 
quantum wells, nominally 39, 30, 23, and 18 ML wide. 
The sample was grown by MBE at 630°C on a two-inch 
undoped GaAs [100] wafer using 8 nm AlAs barriers and 
50 nm GaAs spacers between the wells. This barrier de- 
sign results in a negligible charging of the QWs because 
long-lived unpaired charge carriers can tunnel through 
the narrow AlAs barriers into the GaAs spacers within 
microseconds. Rotation of the substrate was stopped 
only during the growth of the wells in order to achieve 
a continuous variation in well thickness across the wafer 
(along the [110] direction) while maintaining a constant 
barrier width. Growth rates were calibrated using re- 
flectance high-energy electron diffraction on a reference 
wafer. The nominal growth rates were 0.8 ML/s and 0.3 
ML/s for GaAs and AlAs, respectively, and a 30% varia- 
tion in GaAs growth rate was observed across the wafer. 
A standard V/III flux ratio of 8 - 10 was used. At each 
position on the wafer the excitonic resonances of the 4 
QWs of different thickness are spectrally well separated, 
so that by choosing the excitation wavelength wc could 
efficiently select the excitonic response of a single QW. 
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FIG. 9: Spectrally and time-resolved RRS intensity for dif- 
ferent QW thicknesses A as indicated. The spectrally re- 
solved RRS is shown on a constant intensity scale for all A, 
together with the homogeneous linewidth (HWHM) 7ssa(<^) 
given in red. The time-resolved RRS intensity was mea- 
sured by spectral interferometry with 300 fs resolution, and 
is divided through by the exponential decay exp{—2t"/) using 
7 = 25 fieV. Note the break in the time axis. The intensity 
scale is adjusted for each A to the plateau of the data. The 
statistical error due to the finite number of speckles is ±3%. 
The simulated RRS using the ES model (see Section IV. C be- 
low) and the parameters of TableQare superimposed as green 
lines. The intensities have been scaled to match the measure- 
ments, and the time-resolved data have been corrected for an 
average polarization decay rate of 14 fieV. 



The spectrally resolved RRS intensity was measured 
using spectral speckle analysis^, which also determines 
the spectrally resolved homogeneous linewidth 7ssa(^)- 
The samples were kept at temperatures between 1.5 K 
and 50 K. The fundamental hhl-el Is exciton reso- 
nance was excited by optical pulses from a mode-locked 
Ti:Sapphire laser of 0.3 ps pulse length. The diameter 
of the excitation spot was 200 fim. The excited exciton 
density was chosen below 10^ cm~^ to keep the influence 
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FIG. 10: Calculated time- (left) or spectrally (right) resolved 
RRS intensity versus A as indicated. Simulation parameters 
as in tabled] Data are offset for clarity. 



of the density-dependent broadening by exciton-exciton 
interaction weak. The secondary emission was passed 
through an imaging spectrometer and was detected by 
a nitrogen cooled CCD array with a Voigt spectral re- 
sponse of 20 fj,eV Gaussian and 7 fieV Lorcntzian full 
width at half maximum (FWHM). The directional res- 
olution was adjusted to resolve a single speckle, i.e. to 
be better than the diffraction limit of the excited area on 
the samplei^ All presented data were taken through an 
analyzer parallel to the linearly polarized excitation im- 
pinging on the sample in Brewster angle, in a directional 
range centered normal to the sample. The time-rcsolvcd 
RRS was measured using passively stabilized spectral 
interferometrjiSSiSi - a simpler approach as compared to 
actively stabilized versions^^iSS. The time resolution of 
this measurement is 300 fs, while an uncertainty of ±3% 
in the measured intensity is introduced by the average 
over a finite number of speckles. In this case, the diam- 
eter of the excitation spot was 100 /im. 

The measurements were taken at a sample tempera- 
ture of 1.5 K. Different positions on the sample along the 
thickness gradient were investigated. The measured RRS 
for varying fractional ML thickness is given in Fig.|51 As- 
suming a linear thickness variation, the observed lateral 
separation of 3 mm between equivalent fractional ML po- 
sitions of 38 and 39 ML thickness was used to calibrate 
the thickness gradient. The absolute thickness was fixed 
by comparison of the measured RRS spectra with the 
simulated ones at integer ML thickness (see Fig. llOf) . 

The measured RRS shows the formation of ML peaks 



as described in the previous sections. The simulated 
spectra and the dynamics for the parameters of tabled 
are given in Fig. llOl and show a general agreement with 
the measurements. Some discrepancies in the spectral 
shape and in the energy splitting of the ML peaks are 
still present for non-integer ML thicknesses. The simu- 
lated RRS dynamics reproduces the weak amplitude of 
the ML oscillations and their suppression close to inte- 
ger ML thicknesses. The simulated oscillation period is 
about 10 % longer than in the measurements, related 
to the discrepancy in the ML peak splitting. The devi- 
ations could possibly be reduced by further optimizing 
the model parameters. However, one qualitative devia- 
tion between simulation and experiment is obvious: the 
presence of a distribution of polarization decay rates 7 
within the cxciton ensemble. It manifests itself as a non- 
exponential decay of the RRS with time, yielding a pos- 
itive curvature of the RRS dynamics when corrected for 
a single exponential decay. This non- exponential decay 
was already noted earlier—, and attributed to a distribu- 
tion of 7 values within the exciton ensemble. The dis- 
tribution of 7 values is even more directly observable in 
the spectral speckle analysis data, where a spectral de- 
pendence of the line-width is found (red curves in Fig. 

El. 

The distribution of decay rates originates from the in- 
homogeneous nature of the exciton states. Each local- 
ized exciton COM eigenstate is characterized by a differ- 
ent rate of radiative recombination and phonon scatter- 
ing to other exciton states. To model this distribution, 
one needs first to compute the exciton eigenstates - a 
much more computationally demanding task that could 
be avoided in the first part of this work by solving di- 
rectly the dynamical RRS equation versus time, with 
the drawback of having to use a constant value of 7. Once 
the exciton eigenstates are known, the distribution of 7 
can be obtained by computing their radiative recombina- 
tion and acoustic phonon scattering rates in first-order 
perturbationSiiii2£. We have implemented this approach 
for the presented disorder model, as we discuss in the next 
section. The resulting simulations are superimposed to 
the measurements in Fig.|21(green curves), and we discuss 
their comparison with the measurements after introduc- 
ing the model in the next section. 

To evaluate to what extent the present disorder model 
is applicable to excitons in QWs of different thicknesses in 
the same sample - which we assume to feature similar in- 
terface structures - we have measured the photolumines- 
cencc (PL) of the QWs of about 30 and 20 ML thickness 
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TABLE II: Parameters used in the calculations for different 
QW thicknesses. 
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FIG. 11: Measured absorption for 20+ and 30+ ML thick- 
ness. Data are deduced from photoluminescence at 50 K lat- 
tice temperature assuming exciton thermalization to 55 K. 
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FIG. 12: Calculated absorption for 20+ and 30+ ML thick- 
ness. Data are offset for clarity. 



at a temperature of 50 K. For this elevated temperature, 
we can assume a thermal equilibrium in the exciton state 
occupational, which allows us to determine the spectral 
shape of the absorption in a simple wa,y^. The resulting 
data (corrected for the temperature induced band-gap 
shift of 3meV) is shown in Fig. ^2 the absorption, 
up to 3 different ML peaks arc visible. With increas- 
ing ML splitting energy, the ML peaks get better defined 
due to the smaller localization length compared to the 
constant interface correlation length. For 20 ML thick- 
ness, the exciton continuum absorption is merging with 
the upper ML peak, a clear sign that the rigid exciton 
approximationii will not be of good accuracy. The ex- 



citon binding energies Eb for the 30 (20) ML QW are 
12 (15) meV, while the ML sphttings are 4 (10) meV, so 
that already for the 30 ML the disorder potential ampli- 
tude is close to the exciton binding energy. The exciton 
binding energies have been determined by a line-shape 
fit to the absorption of non-GI reference QWa^^ and are 
in agreement with the values reported in Ref.l33i 

The calculated absorption is given in Fig.^] The ex- 
citon Bohr radius ob was scaled relative to the 38 ML 

— 1/2 

simulation according to . All interface parameters 

are as in table ^ except a, which has been scaled propor- 
tional to the ML splitting energy (see Table lll|) . General 
agreement between simulation and experiment is found. 
Quantitative differences are narrower lines and less ML 
peaks in the result of the simulation, compared to experi- 
mental data. The homogeneous broadening 7 at 50 K can 
be estimated to be below 200 /xeV, and should thus not 
be the main origin of the different line-widths in theory 
and experiment. 



IV. MICROSCOPIC TREATMENT OF PHONON 
SCATTERING AND RADIATIVE DECAY 

We now describe the theoretical approach used to sim- 
ulate the RRS signal accounting for a distribution of the 
polarization decay rates. We consider two physical mech- 
anisms. First, the radiative decay of each localized exci- 
ton. Second, the polarization decay due to scattering be- 
tween different exciton states assisted by acoustic phonon 
emission or absorption. The details of how both rates 
arc modeled for exciton COM eigenstates are described 
m Ref.Iil The effect of a radiative renormalization of 
the states that results e.g. in a polarization splittingiSiS^ 
is neglected. The typical splitting energies are similar to 
the radiative decay rates. The splitting can lead to an 
additional slow dynamics of the RRS due to polarization 
beating, mainly relevant on a timcscalc longer than that 
considered in the present analysis. 

In order to compute the rates it is necessary to know 
the eigenstates and energy eigenvalues of the exciton 
COM motion within the disorder potential. These are 
obtained by solving the eigenvalue problem given by the 
COM Schrodingcr equation which reads 

(-|^V2 + y(R)^ ^,(R) = E,^j,{R) , (9) 

where ipjCR) are the COM eigenfunctions and Ej the 
corresponding energy eigenvalues. Eq. (PJ is solved nu- 
merically on the same grid used for simulating the time- 
dependent equation When neglecting decay rates, 
we can express exactly the Green's function of Eq. |^ in 
terms of its eigenvalue decomposition 

G(R,R',i) = ^^;(R')^,(R)cxp (-^i) . (10) 
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From this Green's function, the solution P(R, t) of 
for an arbitrary input field i?in(R, <:) can be computed. 
The advantage of this approach is that we can generalize 
the solution by introducing decay rates 7j for each exci- 
ton COM eigenstate. The corresponding COM Green's 
function can be rewritten as 

G(R,R',i) = ^V;(R')V^j(R)cxp (-^^-7,^) • 
j 

This expression can no longer be related to the solution of 
Eq. However, as we will see later in this section, once 
the decay rates have been accurately modeled, it gives a 
better account of the observed RRS dynamics. As for the 
time-propagation simulations we use a delta-pulse excita- 
tion field at kin = 0, i?i„(R, = EinS{t), corresponding 
to an excitation beam impinging on the sample at nor- 
mal incidence. For this configuration, the RRS amplitude 
resulting from the simulation can be rewritten as 

Kut (k, t) = Yl ^1 (^)^J (^'") '^^P (" ~ ^^*) ' 

(12) 

where 

t/jj (k) = y dRtpj (R) exp(ik • R) , (13) 

is the Fourier transform of the exciton COM eigenfunc- 
tion, A being the simulation area. 

A. Radiative rates 

The radiative polarization decay rates rj axe computed 
from the dipole matrix clement, according to the Fermi 
golden rule. Following Zimmcrmann et ali^ they are 

'>-^i:£^«-^-v2)fe«p^ (14) 

Here A is the area of the simulation domain, k is as usual 
the in-plane wavevector, whereas fcp = y/ebEx/{hc) is the 
wavevector corresponding to the light cone, defining the 
boundary of the in-plane momentum region where an ex- 
citon state can decay into a photon mode radiating in the 
z-direction. Because of energy and momentum conserva- 
tion, excitons at fc > /cq result into a surface polariton 
mode characterized by an evanescent light wave in the 
vertical direction. The quantity Fg is the radiative rate of 
a plane- wave exciton state having zero COM momentum. 
Expression H14(l has been averaged over the azimuthal 
angle and over the two possible polarization directions 
of the emitted electromagnetic field. This average is a 
good approximation, despite the irregularly shaped COM 
eigenfunctions, as the eigenstates relevant to the optical 
response are localized on a spatial range much shorter 
than the optical wavelength in the medium, which for 



GaAs is in the range of 250 nm. Consequently, the wave- 
functions are approximately isotropic for k < ko- For the 
numerical calculations, we have used Fq = 34 ^eV, which 
is estimated^ for a 11 nm GaAs/AlAs QW heavy-hole 
exciton having an exciton binding energy of 11 meV. 

B. Phonon rates 

The phonon scattering rates are computed account- 
ing for the exciton deformation potential interaction with 
longitudinal acoustic phonons. A detailed account of this 
calculation is again found in the Review by Zimmcrmann 
et alm^ and we report here just the essential expressions. 
Fermi golden rule gives the following expression for the 
polarization scattering from state j to state I by emission 
or absorption of a phonon 

27r 

= T ^ ^^'"'^ + " + ^'^^^ ^^^^ 

q 

+nqd{Ei - Ej - hujq)] , 

where q is the phonon momentum in three dimensional 
reciprocal space, is the Bose thermal distribution of 
phonons, ujq = Vgq is the phonon dispersion, Vs being 
the sound velocity. The quantity is the scattering 
matrix element of the deformation potential Hamiltonian 
which, after introducing the factorization ansatz between 
relative and COM exciton motion, is given by 



t?, = ^^^;2^(^i'^^K [DeKA,q.)x{q\\/ve) (16) 

+DhKh{q^)x{q\\/Vh)] ■ 

Here, pm is the mass density of the material. Da are the 
deformation potentials for electron and hole, qz and q| 
are respectively the z and in-plane components of the 
phonon momentum, and rja = (rUe + mh)/ma are mass 
factors. The functions Ka and x s-re overlap integrals 
between the phonon wavefunction and, respectively, the 
confinement functions in the z direction and the relative 
Is exciton wavefunction. The expression 

i^i^iK = J d'R.^ji'^) cxp(-iq|| • R)Vi(R) , (17) 

is the corresponding overlap integral for the in-plane 
COM wavefunctions. After summing Eq. (|16|l over the 
qz variable we find for the scattering rates 

PI, = ^|A£;,,KAi?,OE^^' (18) 

^ q,, «^ 

where AEji = Ei - Ej, qz = ^{AEji/hvs)'^ - q^^ is 

the restriction imposed by energy-momentum conserva- 
tion, Lz is the size of the integration domain in the z 
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FIG. 13: Calculated time- (a) or spectrally (b) resolved RRS 
intensity using the time propagation (black) and the eigen- 
state calculation (red). The eigenstate dynamics has been 
corrected for an average polarization decay rate 7 = 14 ^eV. 
Note the break in the time axis. The spectrally resolved data 
are normalized to the spectrally integrated intensity of the 
38.0 ML data. Disorder parameters are given in table|I]TP. 



direction (which is canceled when plugging in the ex- 
plicit expression for the scattering matrix element), and 
n{E) = {cyL\){E / [k bT) — 1)^^ is the Bose energy distri- 
bution. 

The overall decay rate for a given cxciton state j is 
then given by the sum of its radiative decay rate and all 
out-scattering phonon rates pij 



1] = rj+Pj, 



Pj 



I 



(19) 



where we have denoted the overall phonon outscattering 
rate of the j-th state by pj . For the simulations, we use 
the deformation potentials values Dg = 7 eV, Dh = 3.5 
eV, the mass density value pm = 5.37 g cm~"^, and the 
sound velocity Vg = 5.33 x 10'"' cm s~^, typical for GaAs 
bulk material 1^ The electron and hole effective masses 
were respectively me = 0.07 and ruh = 0.18 times the 
bare electron mass. 



C. Simulations: Ensemble averages 

We can now compare the results of the ES calculation 
(Eq. (HSJ) with the one of the TP calculation (Eq.[T}. Us- 
ing the time-propagation disorder potential parameters 
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FIG. 14: Calculated time- (left) or spectrally (right) resolved 
RRS intensity versus A as indicated. The RRS dynamics has 
been corrected for an average polarization decay rate 7 = 
14 /^eV. Eigenstate simulation with parameters as in Table^] 
Data are offset for clarity. 



given in Tabled we show in Fig.^l the time-resolved 
and spectrally resolved RRS intensity for A = 38.0 
and 38.4 ML calculated with both methods. The time- 
resolved RRS calculated in the TP and ES model are 
nearly identical at early times t < 10 ps. This is be- 
cause all relevant decay rates jj are much smaller than 
the inhomogeneous broadening, so that the nonuniform 
polarization decay among the eigenstates is irrelevant 
(exp(— i(7j — 7)) « 1), and the RRS dynamics is deter- 
mined by the energy-level distribution associated with 
the disordered exciton eigen-states^. For later times, the 
energy-level distribution creates a constant RRS signal, 
given by the incoherent sum of the intensities emitted 
by the individual cxciton states, as can be seen in the 
TP results. In this regime the decay rate distribution is 
relevant since exp(— <(7j —7)) deviates significantly from 
unity, yielding a non-exponential RRS decay. The devia- 
tion from a single exponential decay in the ES calculation 
is similar to the one measured in the time-resolved RRS 
(Fig.O, indicating that the ES model is describing its 
main physical origin. 

Comparing the spectra of the TP and the ES calcu- 
lation, some difference are apparent, especially for ener- 
gies between the ML peaks - the ES calculation tends to 
result in more spectral broadening. This difference re- 
lates to the distribution of phonon and photon related 
dephasing rates. Exciton states energetically in between 
the ML peaks are likely to be excited states of a ML 
trench, which have rather small radiative rates. Since 
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their phonon scattering rates are still smaller than the 
average radiative rates, their dephasing rate is smaller 
than average, resulting in a longer and thus spectrally 
stronger RRS emission as compared to the TP result. 
To reproduce the experimentally measured RRS spectra 
with the ES calculation, we needed therefore to adjust the 
disorder model parameters slightly, as detailed in the ES 
row of tablcQ] The resulting calculated RRS spectra and 
dynamics are given in Fig. 1 141 and show a general agree- 
ment with the measurements. However, the calculated 
average polarization decay rate of 14/zeV is lower than 
the measured one of 25 /xeV. Such a deviation was noticed 
already in previous works^, and can have different ori- 
gins. On the theoretical side, the deformation potentials, 
which enter the phonon scattering rates quadratically, are 
only known up to a factor of two. The radiative rate Tq 
can be influenced by the dielectric surrounding including 
dipole-dipole interaction between the exciton states, and 
also directly by the disorder potential if the rigid exciton 
approximation breaks downii*^ resulting in a modified 
electron-hole overlap. On the experimental side, there 
could be an influence of residual exciton-exciton scatter- 
ing, the influence of spectral diffusioniS, and also a sys- 
tematic error due to the finite spectral resolution. Mea- 
surements of the dephasing time of exciton ensembles in 
similar structures by four-wave mixing^ showed a polar- 
ization decay comparable to the one measured here. We 
have made no attempt to modify the model parameters 
towards a better agreement with the experiment since 
we are mainly concerned with the effect of the disorder 
potential in this work. 

Having at hand the exciton eigenstates with their ra- 
diative emission and phonon-assisted relaxation rates, we 
can investigate in detail their properties at different ener- 
gies in the disorder potential. We do this in the following 
for A = 38.0 and 38.4 ML. The calculated spectral dis- 
tributions of the density of states (DOS) and the phonon 
and radiative rates for low temperature (T=1.5K) are 
given in Fig. 1151 The exciton DOS Dx{lo) is compared in 
the top panel to the normalized spatial potential density 
Dv{uj). For long potential correlation lengths, the lateral 
quantization energies are negligible and Dx{oj) should be 
just given by the convolution of Dy{uj) with the 2D ex- 
citon DOS e{huj - Eo) " 



M 

Tip 



, , M 
Dx{^) « — 



Dv{uj')duj' 



(20) 



This general behavior is visible also in the calculated 
data, but due to the finite lateral quantization energies, 
an additional blucshift of Dy^{uj) is observed. 

The distributions of the scattering rates are more inter- 
esting, as they provide information about the character 
of the states. At a temperature of 1.5 K, the thermal en- 
ergy k-QT K, O.lmeV is much smaller that the ML split- 
ting energy of about 2 meV, so that phonon absorption 
processes are not significant. In this case, the phonon- 
assisted transitions are directed towards lower-lying ex- 
citon states with significant spatial overlap. In a ML 
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FIG. 15: Probability distribution of phonon scattering rates 
Dp{p,uj) and radiative decay rates Dr{r,uj) for exciton states 
in a QW of 38.0 and 38.4 ML thickness. Other parameters as 
in TableE] Logarithmic color scale over three decades. The 
RRS spectra are superimposed as white lines, and the average 
decay rate 7(0;) as red lines. The top plots show the exciton 
density of states Dx{u}) and the potential density Z)v(w). 



island, the local exciton ground state has thus a negligi- 
ble p, while the first excited state acquires some p, and 
the second excited state even more. The radiative cou- 
pling shows the opposite behaviour: The ground state 
has a nodeless wavefunction and thus a large k < kg 
component, resulting in a large r (see Eq. (|14|l '). Ex- 
cited states have additional nodes in the wavefunction, 
and thus suppressed k < components. The calculated 
distributions Dp{p,Lu) and Dr{r,uj) of p and r, respec- 
tively, are shown in the lower four panels of Fig. 1151 For 
A — 38.0, the 39 ML regions have only a small spatial 
coverage, and thus form sub-wavelength sized trenches 
surrounded by 38 ML regions. Most exciton states in the 
spectral region below the 38 ML potential energy i^(38) 
are thus localized in individual trenches of 39 ML thick- 
ness and show the discussed behaviour. The states of 
lowest energy (w < 1.559 eV) are dominantly ground 
states, showing a large r and vanishing p. In the re- 
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gion 1.559 cV< w <1.560eV, the first two excited states 
(p-like) show up, with a 3-5 times reduced r, but a sig- 
nificant p w 3 /icV. For 1.5596cV< oj <1.5605cV, also 
higher excited states are present as indicated by a fea- 
ture of Dp{p,Lu) at p « TfieY, and a broadening of 
Dr{r,uj) of the excited states. The ground states show 
a non-monotonic dependence of r on the transition en- 
ergy: decreasing with increasing energy up to halfway 
to the energy of the 38 ML barrier, and then increasing 
again. This can be understood considering the corre- 
sponding evolution of the wavefunction extension: The 
state energy separation to the 39 ML energy (1.5585 eV) 
can be interpreted as quantization energy Eq of the state 
in the trench. For infinite barriers and a square trench 
of size L, we find Eq = {nh/Ly/M, yielding 0.3 meV 
for L=100nm. Small quantization energies thus corre- 
spond to large trench sizes with large wavefunctions hav- 
ing large k < kg components. With decreasing trench 
size, the wavefunction decreases in size, until -Eq reaches 
about half the barrier height. With further decrease in 
trench size, the increasing penetration of the wavefunc- 
tion into the lateral barrier leads to an increase of the 
wavefunction size and thus of the fc < fcp component. 
This behaviour was previously predicted^2ii£ and exper- 
imentally verified^. 

The energy-dependent averaged lincwidth 

7(^) = {1j)e,=Hu (21) 

is given in the lower panel of Fig. ll5l It shows a depen- 
dence similar to the one measured by spectral speckle 
analysis (see Fig.|3l, with maxima close to the integer 
ML energies due to the large radiative rates of the corre- 
sponding weakly confined exciton states. In this compar- 
ison, we have to keep in mind that the weighting of the 
states in the average lincwidth determined by spectral 
speckle analysis^i^Sii^ is not just uniform, as assumed in 
Eq. , but each state is weighted by its RRS intensity, 
scaling like r'^ /{r+p). This weighting leads to a measured 
dephasing which is typically by some 10% larger that the 
unweighted average 7(w). Yet, this weighting cannot ex- 
plain the observed factor of two between experiment and 
prediction. 

For a direct quantitative comparison of the simulated 
RRS spectra and dynamics with the measurements, we 
show in Fig.ini the ES simulation results for the exact 
ML thicknesses of the measurements superimposed on 
the measured data. The comparison shows good over- 
all quantitative agreement, in the general shape of both 
spectral and time-resolved data and, particularly, in the 
amplitude of the beatings in the time-resolved data at 
non-integer ML thickness. We find the following devia- 
tions significantly above the statistical noise, which we 
believe are not related to a lack of fine-tuning of the 
model parameters: 

i) The oscillation beat frequency and the splitting of 
the monolayer energy are about 10 % smaller in the sim- 
ulation, even though the energy shift due to a thickness 
change of one ML is perfectly reproduced. In the model. 



the smaller splitting for non-integer ML thickness is a 
natural consequence of the lateral quantization in the 
islands of finite size. It is also instructive to observe 
that this effect is less pronounced in the simulation of 
the thinner QWs (see Fig. ll2|l since the in-plane size of 
the ML islands is constant, while the ML splitting in- 
creases as the QW thickness decreases. In the corre- 
sponding experiment instead fFig. lll|l . the effect is more 
pronounced. One possible origin of this deviation could 
be the rigid exciton approximation. Removing this ap- 
proximation, exciton states localized in the small mono- 
layer islands would display a larger electron-hole coulomb 
interactioni2&, leading to a shift of the lower ML peak to 
lower energy. Unfortunately it is numerically prohibitive 
to perform the necessary ensemble averages without this 
approximation. The observed deviation of about 0.3 meV 
is a small correction relative to the exciton binding en- 
ergy of 11 mcV, and lies within the expected magnitude 
of this effect. Another possible source of discrepancy are 
finer details of the heterointerface structure that are not 
accounted for in our disorder model. For example a non- 
Gaussian spatial correlation could give rise to a larger 
spread of island sizes, which reduces the lateral quan- 
tization in the largest islands. Extending the model in 
such direction is a possible matter of future development, 
but we feel it should be supported by more extensively 
sampled experimental data. 

ii) The long-time nonexponcntial RRS dynamics is not 
well reproduced close to integer ML thickness. Gener- 
ally, this dynamics is given by the distribution of de- 
phasing rates across the state ensemble, and also, since 
we excite and detect light linearly polarized along [110], 
by the distribution of the fine-structure splitting and its 
oricntationpiS. We have discussed the dephasing rate dis- 
tribution already in the previous paragraph. The fine- 
structure splitting would lead to an initial decay of the 
signal within a time given by the inverse average splitting 
energy, typically in the 20-50 ps regimaifi. This would re- 
sult in a stronger non- exponential RRS decay, similar to 
the deviation observed. The reduced importance of this 
effect for non-integer ML thickness could be due to the 
alignment of the fine-structure along the [110] direction 
in the (anisotropic) ML islands. A numerical calculation 
of the fine-structure splitting in the present model is only 
feasible in a diagonal approximation in the eigcnstates, 
due to computational limitations. This extension of the 
model, together with a polarization direction dependent 
RRS dynamics measurement, could be an interesting fu- 
ture investigation. 



D. Simulations: Single realizations 

Having described in the previous section ensemble av- 
eraged properties, we turn now to individual excitonic 
states in specific disorder realizations. These results are 
important for comparison with the spectroscopy of in- 
dividual excitonic states performed using high spatial 
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FIG. 16: Excitonic states in a single disorder potential real- 
ization V[(R) of (500 nm)'^ with periodic boundary conditions, 
using the parameters of tableland A = 38.0. Top: Phonon 
scattering (pj) and radiative decay rates (rj) of the exciton 
states j with energy Ej. The numbered states j = 1,2,... 
have ascending energy {Ej < Ej+i). The ensemble averaged 
RRS spectrum is given (blue line) for comparison. Bottom: 
Phonon scattering rates pij from the state j to the state / at 
T=1.5K. Linear grey scale as displayed. 



resolutionMi^Sii&. The simulation area used in our cal- 
culations of (500 nm)^ is comparable to the typical sizes 
of the regions probed in this kind of experiments, so that 
the measured and simulated exciton state ensembles are 
similar. 

We have chosen the thickness A = 38.0 for the simula- 
tions as for this fractional thickness well defined trenches 
of 39ML thickness are formed in the 38 ML plateau. Two 
realizations I & II of the disorder are discussed in the 
following, having trenches of more circular or more elon- 
gated and irregular shape, respectively. The simulated 
potential landscapes V[.ii(R) are shown in the top left 
panel of Fig. El and Fig.^] The remaining panels of 
Fig. El a-iid Fig.EI display the square modulus of some 
exciton COM wavefunctions selected among the lowest 
energy eigenstates. Each of these states is, in general, 
well localized in a local minimum of the disorder po- 



FIG. 17: Disorder potential Vi(R) (top left) on a scale from 
1.5575 (white) to 1.5635 eV (black). a)-e) Probability dis- 
tributions l^&Qp of selected excitonic states as labeled. The 
grey scale is from zero (white) to a) 800, b,c) 500 and d,e) 
300 ^tm"^ (black). 



tcntial. We have therefore chosen to plot in each panel 
several states having similar features, according to the 
discussion below. In the plots, these excitonic states 
are labeled by an eigen-state number j, sorted in order 
of ascending eigen-energy Ej. Each state is character- 
ized by its radiative and phonon scattering rate Vj and 
Pj, which are plotted in Fig.ll6land Fig.EI respectively. 
From the mutual phonon scattering rates pij we can iden- 
tify corresponding ground and excited states of individual 
trenches. 

We start the discussion with realization I. The state 1 
is the global exciton ground state of the system and thus 
also a local ground state. It shows a large r w 10 /icV 
(corresponding to a dipole moment /i w 105 Dcbye as 
given by /i = V dTrhreQC^n^^uj^^ using the refractive in- 
dex n « 3.5), and negligible phonon scattering pi, as 
expected. States 2,3,5,7 show a similar behaviour, and 
should thus be attributed to local ground states. States 
4,6 instead have a reduced r w 3 /ieV and a significant 
p w 3 /icV. They are excited states of state 1, as can be 
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FIG. 19: As Fig. ll9l for the disorder potential realization 
V^ii(R). 



seen from the phonon scattering matrix pij given in the 
lower panel of Fig.^1 which shows sizable pi4 and pig, 
but negligible p24, P34, P26, Pse- From the py, we find 
also that the states 11,14,15 are further excited states 
of state 1, showing also relaxation into the states 4,6. 
Similarly, states 9,10,18,20,23 are excited states of 2, 
and 8,12,17,25,27 are excited states of 3. Energetically 
close to the 38 ML potential the local ground state 34 
is formed, which has an exceptionally large r w 19 ^eV, 
indicating its large spatial extension. 

This discussion can be verified by observing the spatial 
wavefunctions of the various states, which are displayed 
in Fig.lTTI The states 1-4-6 represent a s-px-Py sequence, 
where the s,p,d,e,.. states refer to the sequence of quan- 
tized state groups in a two-dimensional harmonic oscil- 
lator which have the degeneracies 1,2,3,4,.. . The higher 
excited states 11,14 (and 15, not shown) are d-states with 
2 nodes, but are already quite deformed due to the non- 
perfect circular symmetry of the potential trench. Simi- 
lar sequences are observed for the states 2-9-10-18-20-23, 
3-8-12-17-25-27 and 7-16-24-29-38. In the sequence 5-13- 
19-22-31 the asymmetry of the trench is so strong that 
the d and e states are strongly mixed. The state 34 is a 



rather large state localized by a 39 ML trench about the 
size of flB, for which the confinement potential is there- 
fore already reduced due to the averaging over the cxciton 
relative wavefunction |(/)(R')p. 

Generally we can say that for realization I, most of 
the states below the 38 ML peak are well described by 
ground and excited states of single rather round ML is- 
land trenches. In realization II, most of the trenches are 
instead elongated and bending, some even show bifurca- 
tions (see Fig. ll9|l . Accordingly, the states localized in 
the trenches do not resemble the s-p-d state sequence. 
Looking at r and p in Fig.^J the states 1,2,4,6,13,15,16 
are local ground states. State 1 has the excited states 3,5, 
and shares the excited states 8,9,10,11 with state 2. This 
is due to the connected nature of the trenches, which sup- 
port multiple local ground states within a trench struc- 
ture. Such a structure results in a exciton state sys- 
tem with two metastable local ground states (here 1,2), 
and excited states delocalized over both ground states 
(here 8,9,10,11) which allow to coherently couple the 
two local ground states by coulomb and exchange inter- 
action with the excited state. Such coherent coupling 
has been observed in recent literature using nonlinear 
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rate pj for selected exciton states of the disorder potential 
realization l/i(R). 



spectroscop y^Ti^? . 

The different localization of the states results not only 
in different zero temperature dephasing rates, but also in 
a different temperature dependence. Strongly localized 
states show an activated temperature dependence. Espe- 
cially the local ground states (s-like) do show vanishing 
phonon assisted dephasing at zero temperatures, and an 
activated behavior for T > proportional to the phonon 
occupation at the energy of the transition to the first 
activated states (p-like). This is exemplified in Fig.l20l 
for the states l(s),4(p), 14(d) and 34(s) of the first dis- 
order reahzation (see Fig. [T7|l . State 1 shows the largest 
activation energy, given by the distance to its first ex- 
cited states (4,6) of about 0.6 meV. The excited states 
also show an activation on top of the zero-temperature 
offset, but with a smaller activation energy as the higher 
excited states are typically energetically due to the finite 
barrier height. The weakly localized state 34 shows a 
weaker activated behavior, similar to what would be ex- 
pected for an ideal quantum well, where phonon energies 
between about 0.1-1 meV (for a 11 nm GaAs QW) can 
scatter the optically active k = state to exciton states 
of higher /s^SiS. These predictions are in agreement with 
measurements of the temperature-dependent line-widths 
of single exciton states observed in PL experiments^. 
It is interesting to note that not only the activation en- 
ergy, but also the linear slope for higher temperatures 
(fcT > 1 meV) varies between the states. The localiza- 



tion has thus an effect on the phonon assisted dephasing 
even at thermal energies much larger than the localiza- 
tion energies, creating a significant distribution in the 
exciton dephasing rates at elevated temperatures. 
V. CONCLUSIONS 

In conclusion we have presented a model for the dis- 
order of the heterointerfaces in quantum wells that takes 
into account several important aspects of the interface 
formation, like surface diffusion, monolayer island for- 
mation, segregation & intcrdiffusion, and interface corre- 
lation. While not attempting to model the growth pro- 
cess microscopically, the main features observed in struc- 
tural investigations of the heterointerfaces are considered. 
The spectrally and time-resolved resonant Rayleigh scat- 
tering (RRS) spectra calculated for the interface model 
are in remarkable agreement with both the time-resolved 
and spectrally resolved measurements, indicating that 
the main features of the exciton disorder potential are 
captured by the model. Calculating the RRS using ex- 
citon eigenstates and microscopically calculated dephas- 
ing rates due to radiative decay and phonon scattering 
instead of the often used time-propagation model pre- 
dicts the observed spectrally dependent homogeneous 
linewidth and the non-exponential RRS dynamics at 
long times. The calculated localized exciton states in 
the monolayer islands including their phonon scattering 
and their dipole moments are also relevant for compari- 
son with experiments on individual excitonic states per- 
formed using high spatial resolutionMiiSii&. Our result 
indicates that the present disorder model contains all the 
important features of semiconductor heterointerfaces fab- 
ricated by growth-interrupted molecular beam epitaxy, 
and can substantially improve our understanding in a 
large variety of experimental situations. 
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The electron and hole wavefunctions in the in-plane coordi- 
nate R' are oc exp(— M/me,h|R— R'l/as) with the in-plane 
masses m^^h of electrons and holes. For me = rrih, the re- 
lationship used in our work is obtained. In other cases, the 
convolution function should contain a sum of the electron 
and hole probability, weighted with the respective inter- 
face potential (see Ref. llll) . Within the rigid exciton ap- 
proximation, the difference to the chosen approach is not 
relevant. 

In this case, our model results in a vanishing disorder po- 
tential, as the change in the absolute exciton position in 
z-direction is not considered. However, realistically such a 
disorder contribution would exist, giving rise to repulsive 
potentials at the ML steps. We expect these contributions 
to the disorder potential to be weak compared to the ML 
thickness fluctuations as long as their lateral extension is 
larger than the exciton bohr radius. 

As the exciton wavefunction is not constant within the 
well, a weighted averaging should be performed. 



